"0","```r
constrained <- chclust(scpdSM_bray, method=\"coniss\")
bstick <- bstick(constrained, 10)
ngroups <- min(which(bstick$dispersion>bstick$bstick))  # in this case, no clusters should be retained because even with two groups, the broken stick values are higher than the observed Sum of Squares.
x <- strat.plot(scpdSM, yvar=ages, title=\"\", y.rev=T, cex.xlabel=0.5,cex.yaxis=0.5, cex.axis=0.5, clust=constrained) 
#addClustZone(x, constrained, col=\"red\", ngroups)
```"
